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ABSTRACT 



In a previous paper [1], the steady-state behavior of a finite queue 
which accepts batch poisson inputs and received service from servers operating 
in synchronous mode was studied. An analysis was successfully completed via 
the application of the Residue theorem in complex variables. This work 
extends the study in [1] to include the effect of routing and buffer sharing. 
Upon the arrival of a batch each customer determines its route independently 
according to certain probability distribution. Buffer sharing with minimum 
allocation studied in [5] is also considered. Results obtained include state 
probability, blocking probability, delay, and throughput. Validity of 
analysis has been verified by computer simulations. 
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1 . INTRODUCTION 



In a previous paper [1] Chang and Chang study the steady state behavior 
of a queueing system under the following conditions. 

1 ) Customers arrive at the queue in batches according to a Poisson 
process with mean rate A batches/sec. 

2) Each batch carries a random number of customers. The size of each 
batch is a positive integer*‘valued random variable which may follow otherwise 
any arbitrary probability distribution. 

3) The system provides a capacity of accomodating N customers. 

^4) Customers receive services from m servers in the manner of 

f irst>^come*-f irst-served . Each customer requires a constant service time of T 
seconds. The servers operate synchronously in the sense that m customers can 
be removed form the system constantly every T seconds. 

5) Upon the arrival of a batch if the remaining space is not enough to 
accept every customer in the batch then the entire batch will be completely 
rejected . 

The system described above is an approximate model of a packet switch in 
a computer communication network. Although the problem is combinatorially 
very complex, it has been successfully solved via the application of the 
Residue theorem in complex variables. In a separate paper [2], an alternative 
approach called minislot approximation was introduced to solve the problem. 

In a practical system such as computer communication network mentioned 
above, a packet switch may have several output channels and packets may have 
their own preferences in selecting an output channel. In order to make the 
model more practical and the results more valuable, part of the purpose of 
this study is to extend the work in [1] to include the effect of routing. In 
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other words, in addition to the above assumptions, should a batch be 
acceptable each customer of the batch decides from which server to receive 
service independently of the other customers acording to a specific 
probability distribution. We shall use r^ to denote the probability that a 
customer will be routed to server i. Clearly, 0 £ rj _< 1 and r^ + ... + 
rj^ = 1 . This type of routing is usually referred to as random routing in the 
area of computer communications. 

The purpose of buffer sharing among customers routed to different 
outputchannels is to achieve efficient utilization of buffer. Several papers 
C3] ^ [5] have discussed this problem. It is pointed out in [5] that sharing 
with minimum allocation (SMA) performs better than complete sharing (CS) when 
traffic is high. In CS the entire buffer is accesible by all customers. In 
SMA each user has a reserved area which can be accessed only by customers 
routed to that server. In addition to the reserved region there is still a 
shared region. This work also tries to examine the effect of SMA on the 
behavior of a finite queue. For convenience we study the problem in which 
each server only has a reserved space capable of holding only one customer. 
Fig. 1 depicts the conceptual model of such a system. We believe the results 
can be modified to the situation where servers can have reserved areas of 
different sizes. 

The approach using the Residue theorem proposed in [1] will again be used 
to obtain results such as state probability, blocking probability, average 
delay, and system throughput. However this approach has to be modified. The 
modification is necessitated by routing consideration. For example the 
derivation of state transition probability via the Residue theorem not only 
depends on the remaining buffer space but also on the set of busy and idle 
servers. This will be demonstrated in Section 2. 
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Section 3 of this report discusses the extension of this work to unequal*- 
rate servers. Conclusions are made in Section 

Throughout this work ergodicity of the process is assumed so that 
essemble-averages can be replaced by time'^averages . 
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2. MAIN RESULTS 



Let -it [3 denote the probability that the system is at state b = (b^ , 

b^j) at the beginning of a time slot. In b = (b-| b„j) such that b = b^ + 

... + b 0 ), bj^(>_ 0 ) represents the number of customers which are to be routed to 
server i while b denotes the total number of customers in the system. Let x = 
(x-| , ... , Xj„) such that each xj is nonnegative integer and x-| + ... + X[„ > 0 
be an arriving batch of size x-| + ... + x„j in which Xj of them are to be 
routed to server i. Suppose N = 10, m = 2, and b = (1,2). Also suppose 
three arrivals ( 3 , 1 ), ( 2 , 6 ), ( 1 , 1 ) (according to their order of appearance) 
have occurred within a slot. The first arrival (3.1) is acceptable and will 
lead the system state from (1,2) to (^1,3). The second arrival (2,6) must be 
rejected since 2+6+4+3=15> 10. The last arrival (1,1) is acceptable 
since 1 +1 +4+3=9> 10. In other words if b = (1,2), the arrivals 
( 3 , 1 ), ( 2 , 6 ), ( 1 , 1 ) all together result in four and two customers effectively 
accepted by the system and to be routed to server 1 and 2, respectively. We 
use a = (4,2) = (3.1) (1.1) to illustrate such an event and use Pg.^ to 

denote the corresponding conditional probability. In general, a = (ai , ... , 
a^p) in which each a^ is a nonnegative integer such that a = a-| +...+ a[„ > 0 . 

Let n = (n-| , ... , n„,) denote the system state at the beginning of the 
next slot. Then the steady“state behavior of the system is characterized by 



(1) TT = I TTbPa*b 

n (a, b) eAB(n) 

where 



(2) AB(-) = {(a, b) 



0 < (b. 

- 1 



1 ) 



< n^. 
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m 

I (b. ^ 1)^ < N - m, 
i=l 



a + b ^ N, 
m 

and J, (a. + b. *^l)'*'^N‘-m} 

i=l ^ ^ 



In (2), 



( 3 ) 



(x)^ = 



X if X > 0 

0 if X < 0 



Physically, AB(n) represents the set of (a, b) such that the system can 
transit from b to n via the acceptance of a. In (2), (b^ ^ I)"*" represents 
the number of customers (among the initial bj) remain in the system at the 

m 

end of the slot. ^ (b. - 1)"^ must not be greater than N*^m since at most m 

i=l ^ 



m 

customers can be removed from the system per slot. ^ (b. •-1)'*' = N=-m 

i = l ^ 

occurs when b = N, i.e. system is full, at the beginning of a slot and all 
servers are busy during the slot, (a^ + bj^ 1)"^ represents the number of 
customers which stay in the shared region and to be routed to server i. For 
example if bj^ = 0 and a^ = 1 , then this customer will be placed in the 
reserved area and there is no customer waiting for server i in the shared area. 
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m 

Since the shared area is of size N-m, we have ^ (a. + b.^1)'^<N-m 

i = 1 ^ ^ 

in (2). 



in (1) must be solved together with 



(4) 



I 

b e BB 



1 



where 



(5) 



BB = I ^ 



0 < b. < N 
*- 1 ^ 



(m^*1), 1<i<m, 0<b<N 



Physically, BB denotes the set of feasible states. In (5), N ^ (m - 1) 
represents that the system can have at most N “ (m - 1 ) customers to be 
routed to server i, 1 £ i £ m. 

In order to solve (1), we need to find Pg, ^ first. Let Xj^ be a sequence 
of batch arrivals arranged in the order of their appearances in which Xj^ = 

(Xii ,Xi 2 , . . . ,Xim) . We have explained previously that x^j represents the 
number of customers in the ith batch of {x^} which select server j to receive 
their services. Next we use B and I to denote the set of busy and idle 
servers, respectively within a slot. A customer which arrives within a slot 
must wait till the beginnning of the next slot to receive service even if the 
server it selects is idle when it arrives. This is because the servers are 
running synchronously at the same speed. The duration between the arrival of 
a customer and the beginning of the next slot will be called residual period 
of this customer. Let R denote the size of the remaining free buffer 
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including both the shared and reserved region at the beginning of time slot. 

Let I B I and | I | respectively denote the number of busy and idle servers. 
Clearly, R > | I | and R - | I | represents the size of free buffer in the 
shared area. 

The approach using the residue theorem proposed in [1] can be used to 
obtain p^, However it has to be modified. 

Example 1. Suppose R = 5, m = 3, B = { 2 }, and I = { 1 , 3 }. This means 
the shared region can take no more than 3 customers. Thus in {xj^} 

{(i»,3,0), (0,3,1), (1,11,2), (5,1,2), (0,0,1), (0,3,0)} only (0,3,1) is 
acceptable and the rest are rejected. The arrival (0,0,1) is rejected since 
after the acceptance of (0,3,1) not only the shared region becomes full but 
also the buffer reserved for server 3 is occupied. Thus the only customer in 
(0,0,1) which can only be placed in the shared region gets rejected. {(0,3,1)} 
is called an acceptable pattern when R = 5, B = { 2 }, and I = {I, 3 }. {(0,3,1)} 

further divides {xj } into two subsequences Cq = {(4, 3, 0)} and C-| = {(1, 4, 2), 
(5,1,3), (0,0,1), (0,3,0)}. Cq and Ci will be called blocked subsequences of 
{xi {• 

In general, we use {a^ a^j} to denote a general acceptable pattern 

such that = (a. , a, _ . . . , a, ) and a, = a, , + ...+ a, . Let v. and u. 

k k1 , k2, ’km k kl km 11 

respectively denote the number of unoccupied reserved buffers and the number 

of unoccupied shared buffers before the acceptance of a^ + i , i = 0, 1 n^1 . 

For i = n, Un and v^ are defined similarly except after the acceptance of a^. 

Use Rj to denote the size of the remaining free shared and reserved buffers 
between the acceptances of a^^ and + i . Then any arrival pattern {xj} which 
contains al » . . . , an a subset may have the following blocked subsequences 
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Cq t ^1* •••! such th3t ^ ^ ••• ^ijni^ with ^ij * ^ijl 

... + Cj^ jn, satisfies 

(6. a) Coj > Rq = R or the number of unoccupied reserved buffers selected 
by customers specified in £oj is less than Cqj -< Vq for Vq + 1 £ Coj £ 

Oq '^o ~ ^o 



(6.b) 



c. . 
ij 



R. = 
J 



k=1 



or the number of unoccupied reserved buffers 



selected by customers specified in is less than c^j - Vj_ for Vj^ + 1 £ Cj^j 

£ Cjj £ + Vj^ = Rj^, i = 1 n. 

Define for i = 1, ..., n and after the acceptance of a^ 

(7) ®i Ij I server j for which the reserved buffer is occupied by some 
customer } 

(8) li = {j 1 server J idle and its reserved buffer is still empty} 
Clearly, Bj U = (l , 2, . . . , m) , Uj^ = I li | » s^d Vj^ = | j . Next, let Bq = 
B and Iq = I where B and I have been defined previously to be the number of 
busy and idle servers, respectively, at the beginning of a slot. Also define 
Qj (x, Rj^, Bj^, Ij^) to be the probability that an arriving batch of size x is 
rejected because the number of servers in selected by customers in the 
batch is j which is less than x - Vj^. For example 



(9) 



q.(x, R., B., I.) = ( 



Bi 



X > V. + 1 
1 



while 



( 10 ) 



Q 1 ( X , R^ , ® ^ ^ 



2 

k c Z e Bj^ {k} 



' E 

2, G Bj^ 
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In general, 



= ^ C( ®i» Ik})], X > V +2 

r 1 O 11 '^1 

kelj^ 2,eBj^ {k} 



( 11 ) 



Qj (x, R., B., I.) 




c 

kj) eKKj 



< E 

i e BU 





X 



^ t'^1* ••• kJ)]. X ^ V + j + 1 

k=0 ^ 111 j 1 



where 



(12) KK^ = |k^ , 



.j) 



each k„e I . and k, , 
Hi 1 



k. distinct} 
J ' 



Eq. (11) above can be easily encoded into computer program,. Our expreiences 
show the enumeration of qj(x, Rj^, B^, Ij) does not require too much computer 
time . 

Let 6i denote the probability that an arriving batch is rejected between 

the acceptance of aj^ and a^ + i for i = 1, 2 n-1 . 6© and Bn resepctively 

denote the above probability before the acceptance of ai and after the 
acceptance of an- Then from (6. a) and (6,b) 



X - V 



(13) 






X = 



E 

j = 0 



g.( X, 



'O’ °0’ 






- 10 - 



where gx denotes the probability that a batch contains a total of x customers 
and 

(H) i. ■ 8, . , * B, , 2 » ... 

represents the probability that a batch contains more than x customers. 

Next 

(15) 

In general 



X “ V. - 1 

qj ( X, R B . I )] 
j = 0 111 

Let Ng be the random variable representing the number of batches arrive 
in a slot and let AP stand for acceptance pattern then 



( 16 ) 



®i " 



R. 

E* 



= Vi + 1 



g C 

X 



1 X ~ v*! + 1 



X *■ V.J - 1 

®x '• ^=0 '^1’ ®1' ^1^^ 



(17) p[aP = {ai , §2 §n}j Ng = k] 




where a-| + ... + an 1 approach using the Residue theorem [1] we 

obtain 

(18) p[AP = {ai , £2, ... , a^} I Ng = k] 
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= CAP 



1 

27Ti 



c 



dz 



k 

z 



j5o - “j> 



where c > 6o> t ••• Bn 



(19) 



CAP 



n 

n 

i = 1 




m a. . 

n (r.) ] 

j = 1 ^ 



(throughout this report the i assoicated with 2iri denote , otherwise it can 
be used as a subscript or superscript . ) Finally, 

(20) PCAP= 1^1* ^2 * • • • f 1 ^ 



I 

= E PtAP = {ai an! ' = k] 

k=0 



CAP 

2ni 



F(z) 



jio - »J> 



dz 



where = PCN-^ = k] and F(z) = Tq + f-|Z + f2z2 + ... is the probability 
generating function of f|<. The derivation of (20) is similar to that in [1] 
where routing is not considered. 

Define 

(21 ) ya = H< (aia2 . . . a^) 

= { AP = {^1 , . . , = (sji* ®ii> ••• ®im^» 
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t -J. ' s J 

We use I 4'a 1 to denote the number of APs contained in 4'a* 
can be obtained recursively as follows. If AP = (e-| , ... 

(22) AP * en-i- ■] — {^1* •••* e ^ , ®ji+ 1 1 

Let 4* = { AP^, ... , AP|^} where I is an arbitrary positive 
arbitrary acceptable batch, then define 

(23) 4' * g = { AP^ » g, .... APj^ » g} 

The definitions of (22) and (23) are taken directly from [ 



m, 1 5 n < R } 

Similar to [ 1 ] , 4'a^ 
e^} then define 

integer and e an 
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Example 2 



Suppose m = 2 



»(,, 0 ) - 11 ( 1 . 0)11 
’( 0 .,) ■ 11 ( 0 . 1)11 
'^(2,0) " H(2,0)1 , 1(1,0), (1,0)}} 

■ l’(0,0) * '2.0)1 (’'(0,1) * 

'{'( 1 , 1 ) = 11 ( 1 . 1 )} , 1 ( 1 , 0 ), ( 0 , 1 )} , 1 ( 0 , 1 ), 

■ I’(O.O) * '’(I.O) * ">•’>! 

f(2,i) = 11(2,1)}} 11(2,0), (0.1)} , 1(1,0), 

11(1,1), (1,0)} , 1(1,0), (0,1), (1,0)} , 
11(0,1), (2,0)}} 11(1,0), (1,1)}} 

■ '■'(O.O) • '2.’)l “ 1»(2,0) * '“’’’I 

U . (2,0)1 U . (1, 



In general, we have 



( 24 ) 



(a, 



a ) 
m 



( 1 . 0 )}} 

* "•O'll 
( 1 , 0 ), ( 0 , 1 )}} 

1 ( 0 , 1 ), ( 1 , 0 ), ( 1 , 0 )}} 

U . ( 1 , 0)1 

1 )} 
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u 

k w j. '‘v ^2 ^ “2 ■ '<in> * "”* * 

<S '<m>‘ »(a, J„) ' 2 2 ” 



where 



(25) 



(a 



1 ’ 



, a^) - |(k,. 



,M 1 *<i < - 



< k < 
- m - 



but 



(k , . . . , k ) * (0 0 ) } 

I m 

Now let us find P C'l'a]* 

Example 3. Let m = 2. From our previous discussion of (1 3) » (15) ^ (16) 
we have 



(26) 






1 

2t\1 




F(z) 



" “oo> 



dz 



_L [ ^00 ^ 

2 vi J Bqq(z) 



dz 



where Aqq (z) = 1 , Bqq (z) = z - Boo» Sqo obtained from (13) 

assuming Rq, Bq and Iq are known. 

If (a-| , 82 ) = (1,0) we have 



1 / g, F(z) 

(27) f (z - 6 ^^) (z - 6 ,„) 



dz 



00' 



'10' 



— ( 

2 iTi J 



A^q (z) F (z) 

®io 



dz 



1 
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where ^■^Q (z) = giri, (z) = (z " Boo^ “ ^lO^* (27 ), BqO Bio 

be obtained from (13) and (15) under the assumption that Rq, Bq, Iq are known 
and {(1,0)} is only acceptable pattern. 

Similarly if (a-| , a2) = (0,1) we have 



(28) 



’’f’co,,)] - 



1 

2iri 



r Aq^ (z) F (z) 

f “o, 



dz 



where Aqi(z) = gir^, Bqi (z) = (z - 6oo)(z^ “ ^oi ) • In (28 ), Bqo ^nd Bqi 
be obtained from (13) and (15) assuming Rq, Bq, Iq ^ne known and {(0,1)} is 
the only acceptable pattern. If (ai , a2) = (2,0), since 4* (2,0) “ 11(2,0)}} U 
{'^(1, 0) * (1,0)} we have 



(29) P[4' ^2 0)^ = ° 1(2,0)}] +P[AP = {(1,0), (1,0)}] 



1 

2iri 



/ 



A^q (z) F(z) 
B 2 o(z) 



dz 



where 



2 ^20 

A^o (z) - g^r^ ^00^^^ ®00^^^ ■ 



®00^ 



^ 1*^1 ^ 0 ^^^ ^20 
®10 ^ 20 ^ 



(k. 



k2)e 



("^l, 



+ k 



2 ) 



(f) °1< 

(2,0)''l 



+ k. 



A B 

2-k, ,0 20 



(z) 



2-ki ,0 



(z) (z - 6^0 ^ 
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“20 



(z ^ 



«00> “ 



», 0 > 



^20^ 



2 0 

= n n (z « B. .) 
i = 0 j = 0 

Suppose Rq. Bq. and Iq are given, 620 obtained either from 61 of (15) 

by assuming AP = ((2,0}} or from 62 of (16) assuming AP = {(1,0), (1,0)}. 

In general, for m = 2, we have 



f a , a (z) F (z) 

<3<» 2Hf B - (z) 

12 J » ^2 



where 



( 31 ) Aa^ , a^(z) 






rk. + k. 



1’ ”2' " (a , a ) k, 

e 1 2 1 



S "2 

1 ^2 ^(a^^k^, a^-k^)^^^ 



®(a^, a^) 



“(a, »a,a,> 



(32) 







) 

2 



£ BTA (a^ , a^) 



In (32) , BTA (a^ , a^) 





0 ^ ki < a^ , 



0 < k^ ^ a^ 



and all B, 

ki, k2 



must have distinct values}. In other words, there could be several 
combinations of k^ and \<2 so that their values of B are identical. In this 
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case, we select only one of them as representative. This can be explained as 



follows . 

In the process of combining all the related rationale to reach 




(z)/B^ 



(z), it could happen that (k , 
^2 ^ 



k^) • (k^ , k^) but 



I f 

6, . = 6, , . In this case we only let the factor (z - 6, , ) 

k^, k^ k^, k^ k^, k^ appear 



once in B (z). Thus, BTA (a., a^) used in (32) and defined in the line 

a^ , ^ ^ 



following (32) is for this purpose. Notice that . has been defined in 

'^1 * ^ 2 ’^ 

(25). In practice, the value of 6 which appears in B (z) can be 

K^, a^, a^ 

obtained from of (15) by assuming R^, B^, 1^ are given and AP = {(k^, k^)}. 
Finally. 

^(a^ , a^) , b 

" ^ '^'*'(a^ , a^) \ Rq = N “ b, Bq = B, 1^ = I] 

example, all the integrals are assumed to be carrried out 
along a circle whose radius is bigger than any of the poles. 



(33) 



Throughout this 
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I 



(3^<) Pl'Vl = ^ x] 

d ^ a 'I y , djirj J 



-•I 



, k,+...+k k„+...+k k 

■^K + +k( (''') 

2TTi f ®k, • • • m ^ k, ^ ^ k„ ’ "^k ^ 

k <l>a^ 1 1 2 m 



k. k k 

r r . . . r A , „ , (z) F (z) 

^ ^ tn 3^1 krn 

B , 1 / ^ (z “ 6 , 

a. - k , .... a “ k (z) a,a„ ...a ) 

II mm 12 m 



dz 



/ rt 

a.^ , • • . , a 

“b; 

d ^ , • • • ) 



a (z)F(z) 
m 



(z) 



dz 



m 



where K = (k^, ... , k^) and 



<35) a. ■ Z ®K,*... 



1 * • • • r« 

I m 



+k 



k, +...+k k-+...+k k 

( ' ^ ■")...( ") 



k e'f 

— a 



‘1 



m 



k k 

r...r^A , ^,(z)*B (z)] 

' [J] ^1 *** ~ ^ 1 » » » « » 3-m 

- K, a - k “ ®a,, ... a 

11 mm 1 m 



7T ^ ^k k ^ 

(36) B (z) = II “^l’ •••’ m 

‘ * ^m B, , e BTA (a, , . . . , a ) 

k, , . . . k 1 ’ ’ m 

1 m 



In (36) , 
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(37) 



BTA (3>^> •••> ) ** 1^1/ 

I m ' 



k 

m 



0 < k. < a. and all g. 
- 1 - 1 k 



k 

m 



must have distinct values} 



Also, B can be obtained from (15) by assuming R , B , I are known 

K ^ I • • • I U U U 

3nd AP — {(^^> •••t i*6*f 



(38) g. 

^1 ’ 



where 




g,[ E 
J ■ 0 



q^(x, R^, B^, I^)] 



m 

(39a) R, - Ro - E "j 



(39b) = (Bq U {l • u(k^), 2 • u(l< 2 ), .... m • ^(*<jjj) } U {o}} “ {o} 

(39c) I = I {l • u(k ), 2 • u(k-), m • u(k )) 

I u \ id m 

(39d) n 1 I^l 



In (39), u(x) (denotes the unit-step function, i.e., 
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1 if X > 0 

(40) u(x) = 

0 if X < 0 

Notice that |o| and - jo} in B-] is to prevent from including 0 and since 
servers are numbered from 1 to m. 

Finally, p^, can be evaluated as follows 
('»D Pa'H = Bo = B, lo = I. Ro = N - b] 

This is because at the beginning of a slot, the initial state b is completely 
specified by the set of busy servers B, the set of idle servers I, and the 
size of the remaining free buffer space N-b. To evaluate the probability that 
a is accepted is equivalent to evaluate based on B, I and N“b. Once Pa_fti 
is obtained, (1) can be solved. The finiteness of the queue guarantees the 
existence of the probabilities irn regardless of the value of A. The 
minislot approximation porposed in [2] can also be used to obtain Pa.q* 



Now let ir^ denote the probability that the system is at state n at t 



seconds after the beginning of a slot, clearly satisfies 



(42) 



t 

IT 

n 



(a, y e ABT(n) 



subject to 



(43) TT = 1 

n c BB — 

where ABT(n) = {(a, b) | bj £ nj and aj = nj ^ bj| and BB has been defined in 
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(5). Also in (ii2) ^ is the probability | has been accepted into the 

system if at the beginning of the slot the system state was b. Let 0 denote 
the beginning of a slot and define 






* 1 t 

TT = TT I TT 

n ^ J 0 2 



dt 



Then 



(45) 



"6 '’a. b 
(a, b) e ABT(n) 



subject to 



X) 'n ■ ' 

n e BB 



In (45) 



(47) 






dt 



1 

2ni 




A (z)F*'(z) 



dz 



dt 
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1 1 / 

2ni T J B^(z) 






F^(z)dt dz 



1 

2iri 



T 



A (z) 

■yiT 



/-T 



‘'Xt(l-z) 
e at 



dz 



XT 




(z) - 1 ] 



V®’ ^0=^ 



The derivation of (47) is similiar to the derivation of Pa.^. 



We use to denote the probability that a batch arrives at t (0,T) is 

rejected. Also we use B^ to denote the average blocking probability of a 
batch. Clearly, 




x^'vCn)^! 

g [ ^ q.(x, R(n), B(n), I(n))] 



where 



n e BB 



R + 
‘^N-n 



R(n) 

E 

x=v(n)+1 
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(49a) B(n) = {{ 1 • u(ni), 2 • u(n 2 ) m • u(rinj)} U jo}} - jo} 

(49b) 1(n) = {l m} “ B(n) 

(49c) R(n) = N “ b 

(49d) v(n) = R(n) - | I(n) | 

Define a test customer to be a randomly selected customer. Clearly, the 
probability that this test customer is drawn from a batch of size x is xg^/ 



(g^ + 2g^ ••• ). Let B^ be the blocking probability of a test customer t 

seconds after the beginning of a slot and 




then 



(50) 



1 N „ ® 

- 1 E - I E 






n=0 



x=N-n+1 



x=v 



R(n) 
Tn) + 1 






x-v(n)-1 

[ q.(x, R(n), B(n), I(n)][ 

j=0 J 



where 



(51) a = 2 xg^ 

x=1 

represents the average batch size. 

Let Lq denote the average number of customers waiting in the system, then 
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(52) 



^ b E BB a e AA(b)-^ ^ ^ 



1)^ ) 



a,b 



] K 



where AA(b) denotes the set of feasible a if the system was initially at b and 
can be expressed as follows 



(53) 



AA(b) 




m 

0 < a £ N - b, 0 < ^ (a. + b. 

" j=1 J J 



I)"" ^ N r m} 



m 

The reason that ^ (a. + b.-l)‘^,<N-m has been given in the explanation 
j = 1 J J 

of (2). Next let Lq,i^ denote the average queue length associated with server 
k, then 



( 5 ^) 



"q,k 



Z 

b e BB 



[ 



a e AA(k) 



(a, ♦ (b,M)*) 



,bl"t 



We use Xg.k to denote the effective input rate, expressed in number of 
customers per sec, to server k, then 



(55) X = X [Z 8j AR(n, i, k)] 

h€BB - i= 

where AR(n, i, k) denotes the average number of customers accepted by server k 
if the batch contains i customers and the system is at state n when the batch 
arrives. In (55), i ranges form 1 to N - n since it is necessary for a abatch 
not to contain more than N ^ n+ 1 if that batch is to be accepted. AR(n, i. 
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k) can be accepted as follows. First we define an indicator variable IV(k) as 
follows to indicate whether k is in I(n) of (^9.b) 



(56) 



IV(k) = 



! 1 , if k ei(n) 

0 , elsewhere 



Using the notations defined in (M9) we obtain 



(57) 



AR(n, 




y. j(b rj (1 “ r. )^“j, if i - v(-) + IV(k) 
J K K 



v(n) + IV(k) 

E j( ) V 
J.1 J 



(1 ^ 



i - v(n) - 1 IV(k) 

Zi [q (i-j. 

x=0 



R(-), B(-), !(-)]} 



, if v(q) + IV(k) < i < v(n) + 



I(n) 



This can be explained as follows. When i £ v(n) + IV(k) then the batch will 
definitely be accepted. Therefore in this case the average number of 



customers accepted is simply 



I 



j = 1 




r'^ ( 1 - r ) ^ "^ . 
k k 



However when 



v(n) + IV(k) £ i £ v(n) + | I(n) | then the batch of size i could be rejected 
due to improper routing. The term ^ qx(i”j. R(n) , B(n) , I(n)) accounts for 
this effect and must be subtracted. 

Once Lq,j^ g^d Ag,j^ are obtained, the average waiting time of an accepted 
customer routed to server k, can be obtained from Little's formula as 

follows . 



- 26 - 



(58) 



W 



c ,k 




This holds at steady state. 



Then the average delay of an accepted customer routed to server k is 



(59) 




Let Dq denote the average delay of a customer, then 



m 



( 60 ) 



D 



c 



Although the derivation of Dq is done through (58) via Little's formula, 
the average batch delay has to be obtained in a different manner. This is 
of course due to the variation in batch size. But more importantly, customers 
in the same batch can be routed to different servers. This complicates the 
derivation of D^. 

Let ^ denote the average daly of a batch of size 2.. Then 



( 61 ) 



D 



b,1 




where 



(62) AB = l(a, b) • 0 - a + b - N 1, 2^ (a. + b. ^ 1) - N - m} 

J J 
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In (61), the 1 corresponds to 1 slot of service time while the 0.5 represents 
the average residual period. The residual period of a customer is measured 
from the instant of its arrival till the beginnning of the next slot. Since 
arrival could occur anywhere within a slot, the average residual period turns 
out to be 0.5 slots. The KKi in (61) denotes the set of serves which 
guarantees the acceptance of the batch of size 1. Obviously, if there is at 
least one empty buffer in the shared region then KKi = { 1 . 2, . . . , m } . If the 
shared buffer is already full, then KKi is simple the set of servers for which 
their reserved buffers are empty. For 0^,2 we have 



(63) 



b,2 






(a,b) e AE 



b^a,b ^ 



max[a. + (b 

1,2 1 1 



1) + d. + 0.5, 

*^1 



(K^, 1<2)*KK2 



2 2 1 2 



where 6^ denotes the number of customers in the batch of size 2 which select 
server i and 



(64) AB^ = l(a, b) 



0-a + b- N‘-2, I (a. + b.‘^1)'^-N-m} 

j = 1 J J 



(65) KK^ = {(k^ , k^) 



servers k^ and k^ are chosen by the batch and 



m 



V (a +b. + d. -1)‘*’<N^ 
i=1 ^ ^ ^ 



In general, we have 
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+ 



+ 0.5, 



(66) D 



b,H 



E' 

(a,b) e AB 



^a.b 



(k, 



max [a 
k,)e KK^ 



(b, 



1)^ + d, 



where 




1)^ . d 



+ 0.5] r 

*^1 




m 

(67) AB = l(a, b) , 0 < a + b < N - I, V (a + b . - 1 ) ^ ^ N - m| 

t - j=1 J J 



and , 



(68) KKj^ = {(k^, ... , kj^) servers , ... ,V.^ are chosen by the batch and 

m 

E (a. + b. + d. ‘^1)‘^<N*^m} 

111 

1 = 1 

Finally, the average batch delay is given by 



(69) 



1 - B. 



N 

I 

2,=1 






Concerning delay analysis, most papers in the literature concentrate on 
customer delay, mainly for simplicity. For customer delay, relation between 
delay and throughput or input rate can be established usually through Little's 
formula. However, no such relation is guaranteed for batch or group delay. 

The reasons for the problem studied in this report are variation in batch size 
and routing effect. Similar phenomenon is also noticed by Whitt [6] and 
Halfln [7]. 
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The system throughput S is given by 



2 ) ir I 

b e BB k=1 ^ ^ > 

where BB is defined in (5). 

We have carried out extensive numerical calculations based on the 
results obtained here. Fig. 2“5 are part of them. In these examples we 
assume m = 3, N = 8 and ri = r 2 = r 3 = 1/3. Two batches size distributions 
are considered in each of these examples. One is the well known geometric 
distribution with p = 0.5 or a = 2. The other is specified by gi = gg = 

0.267, g2 = g? = 0.133, 83 = 86 = 0.067, and g 4 = gg = 0.033. Since the 
latter has a shape looks like a suspension bridge this distribution will be 
referred as SB for convenience. 

In each figure, the curve labeled by G corresponds to geometric 
distribuion while the one labeled by SB corresponds to the suspension bridge 
distribution mentioned above. In addition to numerical calculations we also 
carry out simulations to support our analysis . We observe in these figures 
that the agreement between analysis and simulation is extrememly good. 

Fig. 2 shows Bg vs X. First, we observe that in both G and SB, Bg is an 
increasing function of X. This is intuitively reasonable. As a matter of 
fact, Bg approaches 1,0 as X approaches infinity. The rate of convergence 

I 

depends on the actual distribution of the size. Second, we observe that SB 
distribution exhibits higher blocking probability than G. This is explainable. i 
Since the size of a batch in SB is restricted to occur in the range between 
one and eight g-| = gg = 0.267 is considerably higher than the rest. On the 
contrary, for geometric distribuion with p = 0.5, g-| = 0.5 and the 
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I 



distribution is strictly decreasing in batch size. This then implies SB had 
higher blocking probability. Third the gap between G and SB closes up as \ 
increases. This is reasonable since as X becomes big enough the system 
occupancy is high and most of the arrivals will be rejected regardless of 
their actual size distribution. 

Fig. 3 and 4 show and respectively, versus X. For G distribution 
we observe both and D 5 increase as X increases. This is reasonable. 
However, for SB, Dq and D 5 first decreases then increases as X increases from 
zero. This can be explained as follows. As X is low, almost every arrival 
can be accepted. Since the batch of size 8 has a good chance of entering the 
system when X is small, the overall delay is high. When X gradually 
increases, only short batches can enter the system thus the average delay of 
the accepted batches drops. As X is high enough, system is full almost all 
the time thus the delay again is high. Notice that the disagreement between 
analysis and simulation is less than 5 %. 

Fig. 5 shows system throughput versus X. Here we observe G has higher 
throughput. This phenonmenon is consistent with Fig. 2. 
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3. Extensions to Unequal Rate Sservers 



Although routing has been considered in Section 2, the servers there are 
still assumed to be operating at the same speed. In practice, servers 
attached to the queue may run at different speeds. In order to move one step 
forward in the modeling of a packet switch, the purpose of this section is to 
extend the results to include servers possible running at different rates. 

In this section we assume when a customer is routed to server k it will 
take the server (or transmitter) seconds to complete its service for this 

customer. Due to the synchronous nature of servers, the behavior of the queue 
can be described on the basis of t‘-second slots where T is the greatest common 

divisor of A-| A^. In addition to the random routing considered above 

we also consider another type of routing called idle^^server-f irst (ISF) 
routing which is believed to offer better performance. In ISF routing, idle 
servers or servers with no customer waiting will be considered first upon the 
arrival of a packet. For both types of routing we have established 
state'^transition equations, obtained state transition probabilities, and 
derived results such as blocking probabilities, delays, and throughputs, etc. 
Based on these results we have also carried out extensive numerical 
calculations. The validity of analysis has also been verified by computer 
simulations. Figures. 6~9 show part of the numerical results. 
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4. Conclusions 



We have in this work studied the effect of routing and buffer sharing 
with minumum allocation on the behavior of a finite queue which receives batch 
Poisson inputs and provides multiple servers running synchronously at the same 
speed. The main contribution of this work is to obtain analytical results for 
system state probability, blocking , delay and throughput. The validity of 
analysis is not only verified by simulation but also supported by intuitive 
reasonings. The work reported here is an extension of [1]. 

We have also extended the model to include unequal rate servers. 
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